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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  new  dynamical  core  that  would  enable  the  Navy 
Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  to  achieve  higher  horizontal  and 
vertical  resolution,  higher  vertical  extent,  increase  meteorological  skill,  and  additional  predictive 
constituents,  all  without  an  appreciable  increase  in  run  time  efficiency. 

OBJECTIVES 

The  objective  of  this  project  is  to  develop  and  transition  semi-Lagrangian/semi-implicit  (SL/SI) 
integration  algorithm  into  the  NOGAPS.  The  SL/SI  numerical  algorithm  is  the  only  numerical 
integration  technique  that  removes  the  CFL  constraint  on  the  system’s  time  step  and  thus  is  the  only 
numerical  technique  that  will  allow  for  increased  resolution,  more  dynamical  variables,  and  higher 
accuracy,  without  a  significant  decrease  in  computational  efficiency. 

APPROACH 

The  work  in  FY  2010  was  to  prepare  the  NOGAPS  SL/SI  dynamical  core  with  full  physics  and 
coupled  to  the  NAVDAS/AR  data  assimilation  system  for  transition  to  the  6.4  Large-Scale 
Atmospheric  Models  Project. 

WORK  COMPLETED 

•  A  two-time-level  NOGAPS  SL/SI  was  developed  and  tested  against  the  three-time-level  NOGAPS 
SL/SI.  While  computationally  more  efficient,  the  results  indicated  that  the  two-time-  level  was  less 
skillful  than  the  three-time-level  version  of  the  SL/SI  code.  Since  the  completion  of  the  two-time- 
level  tests,  several  significant  changes  have  occurred  in  the  NOGAPS  SL/SI  dynamical  core  (see 
next  two  bullets)  and  further  testing  of  the  two-time-level  SL/SI  will  take  place  in  the  6.4  program. 

•  A  further  modification  for  vector  momentum  equation  was  developed,  which  enabled  a  true  time- 
centered  and  therefore  more  consistent  time  discretization  of  the  momentum  equation.  Tests 
indicated  that  this  fix  had  a  positive  effect  on  the  tropical  forecasts  and  reduced  the  120-hour 
tropical  cyclone  track  error  by  80  nautical  miles. 
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including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2Q  jq  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

4.  TITLE  AND  SUBTITLE 

NOGAPS  Semi-Lagrangian 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory,?  Grace  Hopper 

Avenue, Monterey, CA, 93943-5502 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


•  Several  modifications  vertical  interpolation  methods  were  developed  and  tested  in  an  effort  to 
reduce  the  height  forecast  errors  in  the  upper  atmosphere.  These  included  a  correction  to  the 
Hermite  vertical  interpolation  in  the  upper  atmosphere  and  a  non-interpolating  algorithm  in  the 
vertical.  Forecast  tests  over  the  month  of  August  2009  indicated  that  both  methods  reduced  the 
errors  in  the  upper  atmosphere,  but  the  correction  method  to  the  Hermite  vertical  interpolation 
showed  more  overall  forecast  skill. 
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RESULTS 

The  explicit  vector  momentum  equation  in  the  hybrid  coordinate  system  is  given  by 

dV  -  -  - 

- =  -/xF-VO- 

dt 


dP 
c„0  — 
dn 


V  n . 
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where  V  is  the  horizontal  wind  vector,  /  is  the  Coriolis  vector,  O  is  the  geopotential,  n  is  the  terrain 
pressure,  c  is  the  specific  heat,  0  is  the  virtual  potential  temperature,  and  P  is  the  Exner  function: 


P  = 


VPoJ 


(2) 


The  last  term  in  (1)  arises  because  of  the  terrain  following  coordinates  and  is  zero  when  the  vertical 
coordinate  switches  from  sigma  to  pure  pressure  (at  approximately  100  hPa  in  NOGAPS).  The  three- 
time-level,  semi-Lagrangian,  semi-implicit  treatment  equation  (1)  is  given  by 
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where  the  subscript  term  k  indicates  the  vertical  level  and  the  superscript  n  indicates  the  time  level-  the 
(n+1)  term  on  the  computational  grid  and  the  (n)  and  (n-1)  on  the  interpolated  trajectory.  The  two 
bracketed  terms  on  the  right  hand  side  are  treated  in  a  semi-implicit  fashion  -  that  is  these  tenns  are 
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expanded  to  obtain  those  terms  responsible  for  the  fast  gravity  wave  modes.  After  some  algebraic 
manipulations  of  the  definitions  of  the  geopotential  and  Exner  function  (3)  can  be  expressed  as: 


J*v;  -v®;  -c filin’ 

on 

LM 

Y Clkl  ( pVQ’f'  +  [1  -  J3]V@; -  V©; )  -  C;  ( /Nnn+X  +  [I  -  pfjnn~x  -  V7Vn )  (4) 

7=1 

where  The  C’s  are  constant  matrices  and  vectors  and  P  is  any  number  between  0  and  1.  p  -  1/2 
gives  the  classic  centered  difference  approach  to  the  semi-implicit  algorithm,  but  a  stability  analysis  of 
the  non  semi-Lagrangian  Eulerian  advection  in  NOGAPS  indicates  that  this  choice  is  unstable. 

Stability  is  achieved  for  P  >  0.8  .  Tests  of  the  SL/SI  code  indicate  that  this  criterion  is  also  valid  and  so 
in  NOGAPS  SL/SI: 


vT-vV 


2  At 


P  =  0.8  .  (5) 

In  a  pure  three-time-level  semi-Lagrangian,  the  first  three  terms  on  the  right  of  (4)  require  their  own 
separate  spatial  interpolation  (see  Figure  1). 


SL  T rajectory  Calculation 
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Figure  1.  The  three-time-level  interpolation  in  the  semi-Lagrangian  algorithm.  The  (+)  point  lies  on 
the  computational  grid  and  the  (-)  point  is  the  2  At  interpolation  back  in  time,  indication  where  the 
parcel  was  at  t  -  2  At.  The  (0)  point  is  the  location  of  the  trajectory  at  t- At. 
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However,  in  the  original  three-time-level  version  of  NOGAPS  SL/SI  these  terms  are  given  by  a  linear 
interpolation  in  space  of  the  departure  point  (denoted  by  superscript  -)  and  the  arrival  point  (denoted 
by  a  +): 


(7  x  Kf  =  aCf  x  f;t + (i  -  «)(/  x  f;r ,  (6) 

where  a  =  0.8  .  It  was  discovered  that  this  choice  was  not  as  optimal  as  first  appeared  and  lead  to 
significant  errors  in  the  winds.  Figure  (2)  shows  a  comparison  of  the  tropical  cyclone  tracks  forecasts 
for  the  month  of  August  2009  (forecasts  starting  at  0Z  only)  of  the  T3 19L42  NOGAPS  SL/SI  system 
with  a  =  0.8  and  with  a  =  0.5,  together  with  the  non  semi-Lagrangian  results.  The  track  error  is 
substantially  improved  with  a  =  0.5,  which  is  now  used  in  all  subsequent  version  of  NOGAPS  SL/SI. 
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Figure  2.  A  comparison  of  the  tropical  cyclone  track  error  for  NOGAPS  SL/SI  for  2  different  values 
of  alpha,  together  with  the  results  from  the  non  semi-Lagrangian  (Eulerian).  The  resolution  of  the 

forecast  model  is  T319L42. 
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A  persistent  problem  in  the  NOGAPS  SL/SI  has  been  the  large  temperature  and  height  forecast  errors 
in  the  lower  stratosphere  (100  hPa  -  10  hPa).  In  FY  2010  two  approaches  were  tried  to  reduce  this 
error:  (1)  add  an  adjustment  to  the  vertical  interpolation  in  the  lower  stratosphere,  and  (2)  employ  a 
non-interpolating  scheme  in  the  vertical. 

If  we  let  //  be  the  generalized  vertical  coordinate,  then  the  vertical  departure  point  1  is  given  by 

Vd"  1  =  t1a+1  ~  ,  (7) 


where  rf  is  the  interpolated  (Hermite  interpolation)  generalized  vertical  motion  valid  at  time  t.  In  the 
first  approach  for  the  region  above  100  hPa  the  departure  position  is  adjusted  by  a  small  amount  /  : 


where 

f  =  \>f+  if7>° 
[<  /_  if  t)  <  0 


(8) 


(9) 


It  was  found  by  experimentation  that  the  for  the  L42  configuration  of  NOGAPS  that  f+  =  f  =  0.07 

was  an  optimal  fit.  Equation  (8)  is  clearly  an  engineering  fix,  trying  to  account  for  a  systematic  error  of 
the  interpolation  in  the  vertical.  Tests  with  other  interpolating  schemes  without  any  adjustment  (linear, 
cubic,  and  spline)  gave  even  poorer  results  than  the  Hermite,  so  it  was  decided  to  test  this  approach 
with  the  Hermite  interpolation. 

The  second  approach  is  to  rewrite  the  advection  equations,  adding  to  both  sides  a  pseudo  vertical 
motion  that  forces  the  departure  point  (see  Equation  8)  to  be  on  a  grid  point.  So  if  we  write  the 
advection  equation  for  any  field  F  as 


dF_ 

dt 


d„F 

dt 


.  dF 

+  n —  =  R 
dij 


(10) 


where 


dHF 


dt 


represents  the  horizontal  advection  and  R  represents  all  the  remaining  non-advection 


terms,  then  rewrite  the  equation  with  the  vertical  motion  on  the  right  hand  side  and  add  in  to  both  sides 

dF 

a  pseudo  vertical  motion  term  & —  ,  we  obtain: 

dt] 


dF_ 

dt 


d„F 

dt 


dF 

+  ct—  =  R-{i)-ct) 
or/ 


dF_ 

dq 


The  specification  of  the  pseudo  vertical  velocity  is  that  the  departure  point  given  by  is  the 


(11) 
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_  n— 1  M  n+ 1  o  A 

°d  =  t1a  -2o-A? 


(12) 


computational  vertical  grid  point,  closest  to  the  point  given  by  (7).  This  procedure  does  remove  the 
vertical  interpolation  of  the  field  F  ,  but  does  require  an  extra  interpolation  of  the  new  term  on  the 
right  hand  side  of  (1 1).  The  overall  results  does  require  less  overall  interpolations  and  is  more 
computational  efficient  than  the  full  3-d  semi-Lagrangian. 

Figure  3  is  a  time  mean  plot  of  the  1000  hPa  and  the  500  hPa  Northern  Hemisphere  anomaly 
correlation  (AC)  for  T3 19L42  forecast  tests  (August  2009)  of  the  NOGAPS  SL/SI  without  any  vertical 
interpolation  correction  (SL  CONTROL),  NOGAPS  SL/SI  with  the  vertical  interpolation  adjustment 
given  by  Equation  8  (SL  OMT),  NOGAPS  SL/SI  with  non-interpolation  in  the  vertical  (SL  NIV),  and 
the  non  semi-Lagrangian  NOGAPS  (Euler).  While  the  control,  adjusted  interpolation,  and  non  semi- 
Lagrangian  all  show  similar  skill  as  measured  by  the  AC,  the  SL  with  a  non  interpolating  in  the 
vertical  shows  less  skill  throughout  the  troposphere. 
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Figure  3.  The  Northern  Hemisphere  1000  hPa  and  500  hPa  geopotential  height  anomaly 
correlation  for  forecasts  for  the  month  of  August  2009  comparing  the  control  NOGAPS  SL/SI  with 
the  adjusted  vertical  interpolation  scheme  (SL  OMT),  the  non-interpolating  scheme  in  the  vertical 
(SL  NIV)  and  the  non  semi-Lagrangian  NOGAPS  (EULER). 


Figures  4  and  5  are  the  1  hPa  (47.8  km)  height  and  temperature  errors  for  the  4  forecast  test  runs.  It 
should  not  be  too  surprising  that  the  adjusted  vertical  interpolation  scheme  performs  well,  since  it  has 
been  tuned  to  remove  the  cold  bias  seen  in  the  control  NOGAPS  SL/SI.  The  problem  with  the  non 
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interpolation  scheme  appears  to  be  that  it  has  introduced  a  positive  temperature  bias  in  place  of  the  un¬ 
adjusted  cold  bias.  More  work  will  be  pursued  on  different  approaches,  including  increase  vertical 
resolution,  before  a  final  decision  is  made  on  the  best  approach  to  go  forward  into  operations. 
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Figure  4.  The  Northern  Hemisphere  1  hPa  geopotential  height  errors  (mean  and  root  mean  square) 
for  forecasts  for  the  month  of  August  2009  comparing  the  control  NOGAPS  SL/SI  with  the  adjusted 
vertical  interpolation  scheme  (SL  OMT),  the  non-interpolating  scheme  in  the  vertical  (SL  NIV)  and 

the  non  semi-Lagrangian  NOGAPS  (EULER). 
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Figure  5.  The  Northern  Hemisphere  1  hPa  temperature  errors  (mean  and  root  mean  square)  for 
forecasts  for  the  month  of  August  2009  comparing  the  control  NOGAPS  SL/SI  with  the  adjusted 
vertical  interpolation  scheme  (SL  OMT),  the  non-interpolating  scheme  in  the  vertical  (SL  NIV)  and 

the  non  semi-Lagrangian  NOGAPS  (EULER). 

IMPACT/APPLICATIONS 

The  successful  completion  of  this  project  will  enabled  an  significant  increase  in  the  NOGAPS 
horizontal  and  vertical  resolution,  an  increase  in  meteorological  skill  as  measured  by  both  the  standard 
statistical  (NOGAPS  scorecard)  and  synoptic  measures,  and  the  addition  of  more  constituents  (ozone, 
cloud  liquid  water,  dust,  and  aerosols)  without  a  substantial  increase  in  run  time. 

TRANSITIONS 

No  transitions  to  the  operational  code  at  this  point. 

RELATED  PROJECTS 

This  project  is  related  to  the  6.4  Large-Scale  Atmospheric  Models  work-unit,  sponsored  by  PMW-120, 
which  is  responsible  for  transitions  to  the  operational  large-scale  atmospheric  models  run  at  FNMOC. 

PUBLICATIONS 

None 

PRESENTATIONS 

6.2/6. 4NRL/MRY  Review,  March  9,  2010. 
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